Superconducting qubit storage and entanglement with nanomechanical resonators 
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We describe a quantum computational architecture based on integrating nanomechanical res- 
onators with Josephson junction phase qubits, with which we implement single- and multi-qubit 
operations. The nanomechanical resonator is a GHz-frequency, high-quality-factor dilatational res- 
onator, coupled to the Josephson phase through a piezoelectric interaction. This system is anal- 
ogous to one or more few-level atoms (the Josephson qubits) in a tunable electromagnetic cavity 
(the nanomechanical resonator). Our architecture combines the best features of solid-state and 
cavity-QED approaches, and may make possible multi-qubit processing in a scalable, solid-state 
environment. 
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The lack of scalable qubit architectures, with suffi- 
ciently long quantum-coherence lifetimes and a suitably 
controllable entanglement scheme, remains the princi- 
pal roadblock to building a large-scale quantum com- 
puter. Superconducting devices exhibit robust macro- 
scopic quantum behavior 0. Recently, there have been 
exciting demonstrations of long-lived Rabi oscillations in 
current-biased Josephson junctions [2, L|- subsequently 
combined with a two-qubit coupling scheme Q, and in 
parallel, demonstrations of Rabi oscillations and Ramsey 
fringes in a Cooper-pair box [5|, |fj |jj . These accomplish- 
ments have generated significant interest in the poten- 
tial for Josephson-junction-based quantum computation 
Coherence times t v up to 5 /zs have been reported in 
the current-biased devices |2|, with corresponding qual- 
ity factors Q v = T v AE/h of the order of 10 5 , yielding 
sufficient coherence to perform many logical operations. 
Here AE is the qubit energy level spacing. 

In this paper, we describe an architecture in which 
ultrahigh-frequency resonators coherently couple two or 
more current-biased Josephson junctions, where the su- 
perconducting "phase qubits" are formed from the energy 
eigenstates of the junctions. We show that the system is 
analogous to one or more few-level atoms (the Josephson 
junctions) in a tunable electromagnetic cavity (the res- 
onator), except that here we can individually tune the 
energy level spacing of each atom, and control the elec- 
tromagnetic interaction strength. 

Other inv estigato rs have proposed the use of electro- 
magnetic_p, EIlE! El 13 13 El El or superconduct- 
ing [ij], resonators to couple Josephson junctions to- 
gether. The use of nanomechanical resonators to mediate 
multi-qubit operations has not to our knowledge been de- 
scribed previously, although an approach to create entan- 
gled states of a single resonator has been proposed [l9| . 
The use of mechanical as opposed to electromagnetic res- 
onators has the advantage that potentially much higher 
quality factors can be achieved j3|, with significantly 
smaller dimensions, enabling a truly scalable approach. 
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FIG. 1: Left: Equivalent-circuit model for a current-biased 
Josephson junction. A capacitance C and resistance R in 
parallel with an ideal Josephson element with critical current 
Io, all with a bias current It- Right: Metastable potential 
well in the cubic potential limit, showing the barrier height 
AU that separates the metastable states from the continuum. 
Here there are three quasi-bound states |0), |1), and |2), the 
lower two separated in energy by AE. 



Our implementation uses large-area current-biased 
Josephson junctions, with capacitance C and critical cur- 
rent Jo; a circuit model is shown in Fig. ^ The largest 
relevant energy is the Josephson energy Ej = KIo/2e, 
with a charging energy E c = (2e) 2 /2C -C Ej. The dy- 
namics of the Josephson phase difference S is that of 
a particle of mass M = fv^C/Ae 2 moving in an effec- 
tive potential U(S) = —Ej(cosS + s6), for bias current 
s = Ib/Io [2II |22| . For bias currents s < 1, the po- 
tential U(S) has metastable minima, separated from the 
continuum by a barrier AU = £/(<5 m ax) — U(5 m i n ) —> 
(4\/2/3)£j(l - s) 3/2 for s -> 1", as shown in Fig. d 
The curvature U"(S) defines the small-amplitude plasma 
frequenc y lo p = yJ U"(8 m j n )/M = w p0 (l - s 2 ) 1 / 4 , with 
ujpo = \J 2eIo I hC — \/2E c Ej /H. The Hamiltonian for 
the junction phase difference is Hj — P 2 /2M + U(5), 
with P = —ihd/d5 the momentum operator. The junc- 
tion's zero- voltage state corresponds to the phase "par- 
ticle" trapped in one of the metastable minima. 

The lowest two quasi-bound states in a local minimum, 
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FIG. 2: Proposed architecture for two-qubit entanglement. 
The qubits are the two Josephson junctions, coupled to a 
resonator, shown as a disc with a split gate. 



|0) and |1), define the phase qubit. State preparation 
is typically carried out with s just below unity, in the 
range s = 0.95 — 0.99, where U(S) is strongly anharmonic, 
and for which there are only a few quasibound states 
0, . The anharmonicity allows state preparation from 
a classical radiofrequency (rf ) field, as then the frequency 
of the classical field can be set to couple to only the lowest 
two states. In our scheme, by contrast, single quanta are 
exchanged between the junction and the resonator, so 
anharmonicity is not necessary; we find it convenient to 
work with s between 0.5 and 0.9. 

We focus here on coupling a single resonator to two 
Josephson qubits; extensions to larger systems will be 
considered in later work. The two-junction circuit is 
shown in Fig. [21 The disk-shaped element is the nanome- 
chanical resonator, consisting of a single-crystal piezo- 
electric disc sandwiched between two metal plates, and 
the junctions are the crossed boxes on either side of the 
resonator, interrogated by high- impedance circuits Q . 

The phase qubit state |0) of a single junction is pre- 
pared by waiting for any excited component to decay. 
The pure state |1), or a superposition state a\0) + /3|1), 
is prepared by adding a classical rf current J r f to the bias 
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FIG. 3: (a) Phase qubit storage. Solid curves are for | cio (i ) | 2 
and |coi(£)| 2 , dashed curve is analytic RWA for |cio (t) | 2 , dash- 
dotted curve is s(t). (b) Qubit storage with arctangent bias- 
current profile. All other parameters are the same as in (a). 
Solid curves are numerical results for |cio(i)| 2 and |coi (t)| 2 . 



Referring to Fig. [21 the total bias current of junction 
1 is /del + ^rcs, where 7 rcs is the current through the res- 
onator from that junction. A simple model for the res- 
onator allows us to write I les = C Ies (V + h^sbU), where 
C res is the resonator geometric capacitance, the rel- 
evant piezoelectric coupling constant 27], V the rate of 
voltage change, and U the rate of change of the mechan- 
ical strain. The current / rcs is partly due to the capaci- 
tance C res and partly due to the piezoelectrically-couplcd 
strain U. C rcs , in parallel with the junction capacitance 
C, renormalizes the mass M to M = fi 2 C/4e 2 , where 
C = C + C res ■ 

With the resonator coupled to the superconducting 
phase through the voltage V, the Hamiltonian for the 
combined junction-resonator system is H = Hj + H rcs + 



I<t,i(t) = he + cos(u; r ft) + J rf sin(cjrft). Both I dc and SH. Here H les = hwo^a is the Hamiltonian of the iso- 



I T ' ( C vary slowly compared to h/AE. When u> r f is near 
resonance with the level spacing AE/h, the qubit will 
undergo Rabi oscillations, allowing the controlled prepa- 
ration of linear combinations of |0) and |1). 

The nanomechanical resonator is designed with a fun- 
damental thickness resonance frequency cjo/27t ~ 1 — 
10 GHz, with quality factor Q ~ 10 5 - 10 6 . Piezoelec- 
tric dilatational resonators with resonance frequencies in 
this range, and quality factors of 10 3 at room tempera- 
ture, have been fabricated from sputtered A1N |23|, |24| . 
Single-crystal A1N can also be grown by chemical vapor 
deposition |2fij . Our simulations are based on such a 
resonator, with a diameter d — 1.16 fim and thickness 
b = 0.5 fim Such resonators can be used to co- 

herently store a qubit state prepared in a current-biased 
Josephson junction, return it to that junction, or transfer 
it to another junction, as well as entangle two or more 
junctions. These operations are performed by tuning the 
energy level spacing AE into resonance with hioo, gener- 
ating electromechanical Rabi oscillations. 



lated resonator, where we have quantized the resonator 
displacement field with creation (destruction) operators 
(a), and only included the fundamental dilatational 
mode. 5H is the phase-resonator interaction, 



hC rcs bh 33 ■ . t 

SH = — — -6U 5 = ig{a - a 1 )6, 

2e(l-r?) 

where r\ = 0.054 and the coupling constant g is 

_ ft 3 / 2 C rcs h 33 ^ 
9 ~ (1 - n)e y/fmbtP/A 



(1) 



(2) 



For our model resonator g w 0.820 /ieV. 

In the junction eigenstate basis, the junction Hamilto- 
nian is Hj — X3 m em c m c m.) with creation (destruction) 
operators cj^ (c m ) acting on the phase qubit states. The 
interaction Hamiltonian is 



SH = ig 22 (m\5\m') c] n c m -(a - a f ). 



(3) 
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TABLE I: Final state amplitudes Cmnijt/Qa) for phase-qubit 
coupled to nanomechanical resonator. 



probability amplitude 


Cmn 


Im C mn 


| Cmn | 


coo 


0.010 


-0.003 


0.000 


cm 


-0.257 


-0.966 


1.000 


ClO 


0.009 


0.041 


0.002 


Cll 


-0.010 


0.003 


0.000 



TABLE II: Final amplitudes c mim2 „ for state transfer. 

probability amplitude Rec mim2 „ Imc mim2 „ \c mim2 n 

cioo 0.038 -0.013 0.002 

cooi -0.314 0.152 0.121 

coio -0.882 0.422 0.956 
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FIG. 4: Qubit transfer between two junctions. Solid curve 
is |cioo(i)| 2 , dashed-dotted curve is |cooi(£)| 2 , dashed curve is 
I coio (i ) | 2 - Thin solid and dotted curves show si(t) and S2(t). 



res are |mn) = |m)j ® |w) rcs , with energies 
e m + fiwo n, where n is the resonator occupa- 



The eigenstates of the noninteracting Hamiltonian Hq 
Hj + H, 
E 

tion number. An arbitrary state can be expanded as 
= Emn^iWI™) exp(-iE mn t/H). 

The full Hamiltonian is equivalent to a few-level atom 
in an electromagnetic cavity. The cavity "photons" are 
phonons, which interact with the "atoms" (here the 
Josephson junctions) via the piezoelectric effect. This 
analogy allows us to adapt quantum-information proto- 
cols developed for cavity-QED to our architecture. 

We first show that we can coherently transfer a qubit 
state from a junction to a resonator, using the adiabatic 
approximation combined with the rotating-wave approx- 
imation (RWA) of quantum optics [2^]. We assume that 
the bias current s changes slowly on the time scale h/ AE, 
and work at temperature T — 0. The RWA is valid when 
AE and Hujq are close on the scale of Hloq/Q tcs , and when 
the interaction strength g -C AE. At time t = 0, we 
prepare the resonator in the state |0) res . In the RWA, 
neglecting relaxation, we obtain the amplitude evolution 

ihd t c 0n = -iflv^tWe^ci,,.-! 

ihd t c ln = igV^+T (l|<S|0)e~ Mdt c o ,„ + i, (4) 

where lu^ = luq — AE/h is the resonator-qubit detun- 
ing. We integrate to find the reduced density matrices 
Pi{t) (in the qubit subspace) and p T es{t) (in the zero- and 
one-phonon resonator subspace). The junction phase is 
initially prepared in the pure state a|0)j + corre- 



sponding to the reduced density matrix 

H 2 af3* 



PJ(0) = 



(5) 



We allow the junction and resonator to interact on reso- 
nance for a time At — n/ildi where the Rabi frequency is 
fid = (Qq +lu^) 1 / 2 , in terms of the tuned (resonant) value 
£!o = 2.g|(0|<5|l)|/S. After the interval At, the resonator 
is found to be in the same pure state, 



p res (Tr/n d ) 



-a*f3e 



l«|2 



, (6) 



apart from expected phase factors. The phase qubit state 
has been swapped with that of the resonator. The cavity- 
QED analog of this operation has been demonstrated ex- 
perimentally in Ref . |29| . 

To assess the limitations of the RWA, we also numeri- 
cally integrated the exact amplitude equations 



ihc 



m' n' 



'mn\SH\m'n')e i(E " 



,)t/n i 



(7) 



The Josephson junction had parameters corresponding 
Ref. 0, Ej = 43.05 meV and E c = 53.33 neV. We used a 
4th-order Runge-Kutta method with a time step of 10 fs. 
Our main result is shown in Fig. |3| The qubit transfer 
depends sensitively on the shape of the profile s(t), which 
starts at s — 0.50, and is then adiabatically changed to 
the resonant value s — 0.928. We find that the time 
during which s changes should be at least exponentially 
localized. This can be understood by recalling that the 
RWA requires the qubit to be exactly in resonance with 
the resonator (in the Q — > oo limit). Therefore one must 
bring the system into resonance as quickly as possible 
without violating adiabaticity. The power-law tails as- 
sociated with an arctangent function, for example, lead 
to large deviations from the desired behavior, shown in 
Fig. Elb). The result in Fig. I^a) was obtained using 
trapezoidal profiles with a cross-over time of 0.5 ns. All 
quasibound junction states were included in the calcula- 
tion, and convergence with the resonator's Hilbert space 
dimension was obtained. The junction is held in reso- 
nance for half a Rabi period n/Qd, during which energy 
is exchanged at the Rabi frequency. The systems are then 
brought out of resonance. The final state amplitudes are 
given in Table [I] and are quite close to the RWA results. 
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FIG. 5: Preparation of entangled Josephson junctions. The 
solid and dotted lines are si(t) and S2(t), respectively, and 
the dashed curve indicates the probability for the system to 
be found in the state ( 1 100) - |010))/a/2- 

To pass a qubit state a|0) + from junction 1 to 
junction 2, the state is loaded into the first junction and 
the bias current changed to bring the junction into res- 
onance with the resonator for half a Rabi period. This 
writes the state a\0) + (3\1) into the resonator. After 
the first junction is taken out of resonance, the second 
junction is brought into resonance for half a Rabi period, 
passing the state to the second junction. We have simu- 
lated this operation numerically, assuming two identical 
junctions coupled to the resonator described above. The 
results are shown in Fig. 0]and Table ITT1 where c mim2n 
is the probability amplitude (in the interaction represen- 
tation) to find the system in the state \rn\miri), with mi 
and 7712 labelling the states of the two junctions. 

We can prepare an entangled state of two junctions by 
bringing the first junction into resonance with the res- 
onator for 1 / 4th of a Rabi period [3(| , which, according to 
our RWA analysis, produces the state (| 100) - |001))/\/2. 
After bringing the second junction into resonance for half 
a Rabi period, the state of the resonator and second 
junction are swapped, leaving the system in the state 
(|100) - |010))/V5 with a probability of 0.987, where the 
resonator is in the ground state and the junctions entan- 
gled, as demonstrated in Fig. [SJ Using the cavity-QED 
analogy, it will be possible to transfer the methodology 
developed for the standard two-qubit operations, in par- 
ticular controlled-NOT logic, to this system, using mostly 
existing technology and demonstrated techniques. 
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